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Abstract 

The methodology end . rigorous solution fdnnnWion » pre^eW for ^ 
intensity factors (SIP's, It) end total strain energy release nta >(®BK. of. 
multicracked plate, that has fully interacting cracks and is subjected to a for-fi 
arbitrary stress state. The fundamental perturbation problem is derived, and 
steps needed to formulate the system of singular integral equations whose soluti 

^ to the evaluation of tiTsiF’s are identified. Parametric studies are co^ 
ducted for two, three and four crack problems. The sensitivity and e 1 

of the model is demonstrated. 


1 Nomenclature 
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- inclination angle between inner tips of two parallel cracks 

- direction cosines between two local coordinate systems 

- strain tensor 

- four roots of the characteristic equation 

- angle defining orientation of local coordinate system 

- far-field and total stress field, respectively 

_ components of stress in global coordinate system ^ 

- stress from the perturbation problem in p* local frame 
_ transformed crack stressed to the p* local frame 

- normalized real variables 

- Fourier transform of the stress function with respect to x 

- angle between L — T and X — Y coordinate systems 


1 



a,b,c,d 


bn, 12 >22 >16 >26 >66 

fni 
ki,k2 
ke r 7 

Pi 

9i 

H 

rjx> r jY 

s 

t 

u, V 
w 

(*,,»), (XX), (l,T) 

w 

Bit >12 >22 ,66 

C m 

Dj, Sj, Rj, Qj 

Dh > 

Bll, Ett> Glt, vlt 
Fj(xi,Si) 

a T 

K.„ 

m 


— roots of the characteristic equation (real numbers) 

— half crack length 

— co mpliant matrix coefficient in Xj — Vj frame 

— auxiliary functions 

— mode-I and mode-II stress intensity factors 

— FYedholm kernels 

— normal traction at crack surface 

— shear traction at crack surface 

— position vector defining the origin of a local coordinate system 

— components of the position vector ry 

— Fourier variable 

— real variable along a crack axis 

— displacement associated with x and y coordinates, respectively 

— weight function . . . 

— local, global and material coordinates (L - strong direction) 

— kernel matrix 

— compliant matrix coefficient in. L — T frame 

— functions of s in Fourier space (i.e., constants in x, y-real space) 

— constants of substitution 

— horizontal and vertical distances between crack tips 

— material’s parameters in. L — T frame 

— modified stiffness parameters 

— Airy stress function 

— total strain energy release rate 

— discrete auxiliary function 

— effective stress intensity factor 

— loading vector 

— (iso) stands for isotropic singular, (reg) for nonsingular part 


2 INTRODUCTION 

Consider multiple cracks embedded in an infinite anisotropic plate (Fig. 1(a)). The 
plate is under a far-field stress denoted by o° jX , (in particular ^xx, fxY'™** 

(X Y) is the global coordinate system), and the cracks are defined in their local frames 
( x , y ,) (Fig. 1(b)). The origin of each local frame is defined by the position vector,, 
and the orientation of the local frame with respect to the global frame is defined by the 
angle Each crack is symmetrically situated within its own coordinate system and is 

2 a,- long, as shown in Fig. 1(b). ... a . . 

The general solution formulation can be outlined in four basic steps. e st ^ 
is to derive the local stress equations for each crack in its respective local coor a 
system. This derivation is achieved by defining the fundamental problem; that is a single 
Jack in an infinite anisotropic plane (Fig. 1(b)). The fundamentid problem is then 
decomposed into two subproblems: the problem of the undamaged plate containing an 
imaging crack (Fig. 1(c)), and the perturbation problem (Fig. 1(d)) of a plate vuth 
a single crack subjected to the appropriate crack-surface tractions which are found fro 
the solution of the complementary undamaged problem. The analysis of the perturbatio 
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problem leads to singular stresses that govern local crack-tip behavior. 

The second step is to formulate the total perturbation stress field for each crack, which 
includes the interaction of all cracks through the summation of the transformed lo 
stresses of all other cracks. In the third step of the formulation, the total stress equations 
are normalized. A set of Cauchy-type singular integral equations, expressed m terms 
of unknown auxiliary functions, is obtained by subjecting the total perturbation stress 
equations to the crack-surface traction field at each crack location. The fourth and 
step of the formulation is to express the stress intensity factors (SIF’s, krand k 2 ) in terms 
of the discrete auxiliary functions H vj (r p ) evaluated at each crack tip. These discrete 
auxiliary functions are obtained by implementing of the Lobatto-Chebyshev collocation 
technique. Finally, the strain energy release rate (SEER, G r ) is calculated in terms of 

SIF’s. 



Figure 1: Muliticracked plate geometry and method of solution, (a) Multicracked plate, 
(b) Fundamental problem for jth crack, (c) Undamaged plate, (d) Perturbation problem. 
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3 LOCAL STRESS FORMULATION 

Consider the fundamental problem (Fig. 1(b)), which is deBned as a « «■* ® 

an infinite anisotropic plate; its solution can be obtained by dKompoeing^ 

u- d problem (Fig. 1(c)) and a perturbation problem (Fig. 1(d)). The essence 

of this decomposition is that the traction forces applied ulong ' th ? 
perturbation problem are the oppodte of the obtained stress field of U> * 
at the partic ular location of the imaginary crack. As a result, the undamaged p 
traction field can be defined in terms of the normal fa) and shear fa) stress componen s 

along the imaginary crack surface: 


where 


( 1 ) 
( 2 ) 

(3) 

(4) 

The mixed boundary conditions for the perturbation part of the fundamental problem 
(Fig. 1(d)) are expressed in terms of stresses 

(5) 


Pj( x j) — °i w( x i>®) 

Qj( x j) = (T x j vj{ x j>®) 

^(^,0) = a^ x sinVi + ^r cos Vi-^ sin *Pi 
<r XjVi (xj,Q) = ~ — sin 2 yj + a° XY cos 2 <pj 


J ViVi 


* ~Pj(Xj) and 0‘s ijif — ~ < b( x i) 


along the crack surface (in, * = 0 ond - a, < 4 < <*), -d in tarms of continuity of 
displacements md (6) 


outside of the crack (i.e., * = 0 and |*,| > «„ see Fig. l(b)). Here 

indicates the value of displacement at a point approved from the P«rtm snde i of 

plate, (i.e., y > 0), whereas indicates the same point approached from the negative 

^atili°L the preceding two-dimensional anisotropic plate problem 
can be expressed in terms of the Airy stress function *)(*;, Vi) as 

(7) 


&F , #F . & 1 L + ~!?L = q 

+ 71 da?dy + 72 dx?dy* ^dxdy* W 


where 


7i 


. 

“ » 


72 


= k and 74 = fe'. 


( 8 ) 


Here for general two-dimensional anisotropy &n,&22,&i2,&i«,&26,and *** •» independrat. 

InZt paper, we will be primarily concerned with applications mvol^ umdi^- 
tional fiber reinforced composites; these can be idealized at the macrolevel as a pseudo- 
homogeneous transversely isotropic material. Consequently, the preceding ^ inde- 
pendent constants now become dependent on the four independent elastic chants, 
E ll , Err, G l t and v L t , where (L,T) is the material coordinate system rotated by the 

angle fi with respect to (X,Y), see Fig. 1(b). 
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Thus, 


> 

Bl1 cos 4 (ft - v>) + (2Bi2 + Bee) sin 2 (ft - ?) co s2 (° “ ¥>) + ^ 22 ^ 

2" £• K - 8 + U; + JW -’(« - d«*g» - *>) + Busman - V) 

B12 + (Bn + B22- 2B12 - Bee) sin’(ft - <p) cos (ft - v) 

ySXfo W(a J $ + «"■*! ' *> 

[cos 2 (ft - <p) “ Sin 2 (ft - <p) - § (2B W + Bee) cos 2(ft y>)] sin ( <p) 


(9) 


611 = 

622 = 

5 W = 
bee = 

6ie — 

626 = 

wherc ft.-* «-■* “ d *“* fc Jitv 

Note that this special case of transverse isotropy does not torecover 

l^^^ontT^erimSy iTen^fythe independent constants used in 

^rigorous solution for this stress fiinction can be obtained by employing the Fourier 
transform. Assume the stress function to be expressed ss 

( 10 ) 


^)-s££ <w ~ v ~* 


Then, on substitution into eq. (7) the characteristic equation is obtained. It has four 
complex roots, which take the following form: 


fix = a + ify /X2 — c + td; 
= -a + it>; |i 4 * -c + id; 


(U) 


Wh Tlie A^'stiei function must also satisfy the physical requirrauent that the stress 

function tefbute throughout the domain of the plate. Wore, the 

which are automatically bounded at infinity, can be used for the upper half 

plane (for y > 0), 

F(li #+) 1 jf~ [ Cie («-W)» + e-^ds (12) 

and for the lower half plane (for y < 0), 

F{Xf y -) = ±-J~ [GP*** + C A e^ + ^] e~ Ux ds (13) 

XT A r- f fe r fnr i = 1 2 3 and 4 are functions of the Fourier variable s and are 

determined by using the local stress continuity con^^^ ^ 

half planes ( y = 0) and by using the perturbation boundary 

the are calculated by usingthe second 

denvati^ofThrless Actions [1). Wore, the stresses for the upper half plane are 

M = _L r [c, (-o|s| + its) 1 fMtr + ft (-c|s| + ids) 1 #»-**•] e-^d, (14) 
** 27 r 7-00 L 
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0 .(+) L f°° 5 2 [ + Ctff*** e “* ds (15) 

w 2 TV J-ao *- 

„( + ) = ± f" j»fc 1 (-a|«|+»»)e‘^^ + G(-c|*l+««)' ( “ r ' el * l) ’]«' i " <fa < 16 > 

** 27T 7—oo *■ 

and for the lower half plane 

„£> = -L r [Oj (a|s| + <M J e (* 64+Ql *'J v + C 4 (c|»| + e-“ci> (17) 

a (-) = _J_ /‘°° s 2 [C 3 e (iia+0,,1), ' + C 4 e** d ' +C '*^] e'^ds (18) 

w 27T 7 -oo ^ 

„(-) = i- r « fG, (a| >| + ibs) P*4*» + C. (c|»| + ids) e'""' 1 * 0 *] e'^is (19) 

** 27T 7-00 

The continuity conditions for local stresses <r w and <r„ are identically satisfied, given 

Cj + C^Cs + Ci ( 2 °) 


and 


Ci (bs + ta|s|) + Ci (ds + ic|s|) = C 3 (bs - ia|s|) + C 4 (ds - tc|s|) (21) 


respectively. 

The solution of equations (20) and (21) for C 3 and C A in 
written in the following form: 


terms of Ci and C% can be 


Ci = S\C\ + S 2 C 2 (22) 

C A — S 3 C 1 + S A Ci 


where 


_ |a|(a + c) + i(d — 5)a 

1 ~~ |s|(c - a) + i(d — b)s 
2c[sj 

$ 2 ~~ |s|(c — a) 4- i(d — b)s 
2a]sj 

03 = 


(23) 


S 4 = - 


\ s 

(c-a) 

|$ 

(a + c) 

M 

(c-a) 


-b)s 
— b)s 


The strains are calculated by using the generalized Hooke’s law. The normal strains 


ElX — b\\<Txx + bllPyy + blB&xy (24) 

£yy — biiffxx + b22 a w + btt (T *V 

Rom eqs. (24) the strains for the upper and lower half planes can be obtained. Then 
by using the strain-displacement relations [1]; the displacements for the upper and lower 
half plane can be obtained: 

u(x,y) = / £xxdx ( 25 ) 

v(x,y) = J £ n dy 

To obtain the singular integral equations we introduce the following auxiliary func- 
tions: 
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/,(*)= [u + (*,0)-u-M)] (26) 

Expressions for the unknown constant C, end C 2 can be determined in term, of te 
amdluJry functions, since we know that AC*) and A(*) nonzero only w,thm the crack 
region (i.e., -a < t < a). Therefore, 

f D,fM - DiMt) (28) 

1 J—a D\D^ — DiDz 

_ r* £> 3/1 (*) - (29) 

C/2 " J-a DiD 4 -D 2 D 3 

w ^ ere Di = 2a|s| [|s| (a + c) + i (d — 6) s] bn 

D 2 = 2c|s| [Vl (a + c) - i(d - b)s]bn 

UU \s(d-b)-j(a + c) Ml (30) 

03 = 2a (a 2 + 6 2U c | s | + ids) 

_ lsUf-3(d-6)-»(a + c)M s 

04 - 2c ( c 2 + d 2 )(a|s| + t6s) 

Similarly, C 3 and C 4 can be expressed in terms of the auxiliary functions by using eqs. 

^“e^lTns'fS £ constants C m into the local stress ^(H) 
through (19) results in the formulation of a set ofdouble ^th respect 

to the Fourier variables s (-00 < s < oo) and t ( a -)• ^ £ or 

to s will give a set of singular integral equations with respect to t, which 

any jth crack within its own local coordinate system (x j} yj) : 


a ?m = S LI, [ /ilfe) 6^ + 

a %, = h + 


&22Q0 

Qa 

b&Qo 


where 


Q 0 = ac f(a + c) 2 + (6 - d) 2 ] [(*i - *i) 2 + “ x i)% + + ^ 

x [(tj - XjY + 2d(t j - Xj)yj + (c 2 + d 2 ) j£] 

Qi = i?! (tj - Xjf + 2/i [^2 fa - Xjf + VjR* (*> ~ x i) + K?** 4 ] 
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R 1 = a 2 bc + 4- 2ab<? + 2a 2 cd + <u?d + ad 3 

R t = a 4 c + 2a 2 6 2 c + 6 4 c + 2a 3 c 2 + 2a6 2 c 2 + 2a 2 c 3 + ac 4 

+2a 2 6cd + 2b*cd + 2abc 2 d + 2 a 2 ctf + 2 ac*d 2 + 20 W 3 + ad 4 
i? s = a 2 b<? + 6 3 c® + 2a6c 4 + 2a 4 cd 4 Aa 2 b 2 cd 4 2b A cd 4 a 3 c*d 

+a fc 2 c 2 d + a 2 bcd 2 4 b 3 cd 2 4 4a&c*d 2 4 aV + a^d 3 + 20 M 4 
H* = (a 2 + 6 2 ) (c 2 4 d 2 ) (o 2 c + b 2 c + ac 2 4 ad 2 ) 

Q 2 = Hs [fo - xj) 3 i?e + yi [i?r (*; - Xj? ~ VjR* (** ~ x i) + } ( 36 > 

-with 

Rs = (a 2 + 6 2 )(c 2 + d 2 ) 

He = (a 2 c 4 6 2 c 4 ac 2 4 ad 2 ) 

H 7 = a 2 bc 4 b 3 c 4 2 b 2 cd + ac 2 d + 2aW 2 4 ad 3 

H* = a 3 c 2 + a!> 2 c 2 + a 2 c 3 -6 2 c 3 -2a 2 6cd-26 3 cd 

— 2a6c 2 d - a 3 d 2 - ai 2 d 2 + a 2 cd 2 - b 2 cd 2 - 20 W 3 
H 9 = (a 2 + b 2 )(<? + d 2 )(bc + ad) 

Qs = -Rio(ti - *j) 3 + Vi [“ ( f i “ x i ) 2 + Vi (*i ~ x i) *“ + ^^ 13 ] 

with 

i?io = 6c + ad 

H n = (a 2 -6 2 )c-26d(a + c)+a(c 2 -d 2 ) 

H 12 = 6c 3 + a 3 d 4 ab 2 d 4 bca 9 4 26d (6c + ad) 

«i 3 = (a 2 + ft 2 ) (c 2 + d 2 ) (a 4 c) 

Qa= R* {(tj - xjf R u + Vi [*« (*i - *i? + yiiPi ~ x i) R * + ^ Hl7 ] } ( 38 ^ 

with 


i?i4 = a 4 c 

His = 2ab + bc + ad + 2cd 

r 16 = a 3 + ab 2 + 2a 2 c + 2ac 2 + c* + 2aM + 2bcd + cd 2 

Hit = 2a6c 2 4 6c 3 4 a 3 d 4 at 2 d 4 2a 2 cd 4 feed 2 

Qs= Risitj-Xjf + yjlRnitj-xtf-yjitj-xJRto + yfal (39) 
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c (a 2 + &*) + a (c 2 + rf 2 ) 
a 2 6c + 6 3 c + 2& 2 cd + ac 2 d + 2a6d 2 + ad 
o 3 c 2 + ab 2 c 2 + a 2 c 3 - b 2 c 3 - 2a 2 bcd 
—2b z cd - 2dbc 2 d - a 3 d 2 - aft 2 d 2 + a 2 crf 2 - 6 2 cd 2 - 2 aM 3 


g 6 = ils {(«,- - x,) 3 i?2i + yj [Rn to ~ + Vi to ~ *** + ^ i?13 l } (42) 


R 21 = ad + be 

r ^ = - a 2 c + b 2 c — ac 1 + ad 2 + 2abd + 2bcd 

= fee 3 + a 3 d + ob 3 d + 2i?cd + 211 W 2 + bed 2 

Note, that the special case of an isotropic material can be recovered by the following 
substitutions: b = d - 0, and a = c = 1, thus giving «! = * = * “ * = *“ “ ™ I 

= -Rm = 2; and Rj = Rie = 6. The parameters Q, then become 

Qf’= 4 [(»» -*»)’+«?]’ (43) 

Q<<~ ) = 2% [3 (t/ — ®y) 3 + S$] (44) 

Qi*~> — 2(t J -i i )[fe-%) I -Sj] (45) 

Qf’ % [(t, - Xi? - yj] (46) 

Q< < “>= 2(t i -lj)[(t i -il) 1 + 3^] (47) 

<?'“•' = 2(t J -n J )[(ti-*y) 2 -S?]=<3f“ > (48) 

<?<<“>= -%[(‘i-i) , -^]=0 3 ‘“ > <49) 

and the stresses reduce to the isotropic stress formulas derived previously in reference 

,2 ' This completes the formulation of the fundamental problem (or local stress state) 
for the jth crack. Henceforth, the formulation of the multiple crack problem 
addressed. 
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4 TOTAL STRESS FORMULATION 


The total stress state for the pth crack is defined as the local stress state ofthe 

pth crack (o* x ) plus the contribution to that stress state of all remaining cracks. This 
may be represented mathematically as 


n-1 

po^ z (xp,y p ) = 0?z(*p>i/p) + 53 Hj[ x j( x p>yp)>yj( x p>yp)] 

i=l 

fo r p = 1 , ... , n, where a 3 n is defined through standard tensor transformation of the stresses, 
i.e., 

&rz — PlrPmzGlm 

and # r , Pm* are the direction cosines between the (x j} y,) and (a* ,%>) coordinates with j 
identifying the remaining cracks. Note, this statement does not imply that the concept 
of superposition has been invoked, since the stress perturbation boundary conditions (see 
eqs. ( 5 )) have not yet been utilized to determine the unknown auxiliary functions. 

For functional compatibility within eq. (50), coordinate transformations must be 
simultaneously applied to all remaining jth crack coordinate variables. As a result, the 
dominant part (i.e., the first term of eq. (50)) possesses a singularity whereas the regular 
terms within the summation lose their original singularities and yet still contribute to 

the total stress state, as one might expect. 

By replacing subscript j by p in eqs. (33 and 32), respectively, and evaluating them 
at 3 fe = 0 , the singular terms of the singular integral equations are obtained from the first 
term of eq. (50) applied for shear (o**) and for normal stress (<r w ). Finally, the variables 
x and t are normalized by using and ^ = OpT , where £ and r are defined between 

-1 and 1. Therefore, 

*-&**&*■*&&* " 
where E~ and E w are material-related coefficients proportional to local x-direction stiff- 
ness (denoted by superscript 1 ) or local y-direction stiffness (denoted by superscript 
2) with respect to the local crack coordinate system. Consequently, we will call them 
modified stiffness parameters (MSP’s). The MSP’s are 

™(i) |c(a i + t 2 ) + «(c*+<* 2 )l ( 53 ) 

■» 2 ac [(a + c ) 2 + (6 — <f) s ] &11 


(g+j fl(c» + d»)H + tc) 

- 2ac[(a^cf + (b-df\bn 

™(i) (bc + ad) — 

m 2oc[(o + c) s + (t-<*) J ]6ii 

2) _ (a» + 5?)(c»+rf»)(a + c) 

” 2ac [(a + cf + Ci-d) 1 ]^ 


(54) 

(55) 

(56) 
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Figure 2: Geometric relationships between a pair of cracks and their local variables. 


For the isotropic case 


E* 2) = e$ ) =q 


£P(i) = £p<2) _ tL (57) 

whereas in the orthotropic case, only 

£*( 2 ) = E g) = Q 

The regular terms of the singular integral equations are obtained by transforming 
the renmining stresses into the local pth crack coordinate system sinmltmeously mth 
coordinate transformation. The coordinate transformation between the xj.yj and x r ,y r 
systems is determined from the following geometric relationship (see Fig. 2): 

Tjx + Xjcostfj - yjsintpj = r pX + x p cos<p p — y p sirup p (58) 

rjy + XjSintpj + yjcosipj = r p y + x p sirvp p + y v cos<p p (59) 

where r iX , r jY are the rectangular components of the jth crack position vector i referred 
to the global coordinate system X-Y, and <pj is the angle of rotation between the global 

^On^omp^nt of the regular part of the total stress is obtained by tr ^sfoimation 
of the stresses from jth crack local coordinate system into the pth crack local coordinate 
system simultaneously with the coordinate transformation and substituting y P - U. 

•U - - « - •&)"» *i +a *»d c °* e - si f^ «»> 

= ai 3) r . sin 2 9 + <J^L cos 2 6 - 2a^ Vj sm 6 cos 9 
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( 61 ) 


where 6 = (fi p — <pj. Therefore, we can obtain from eqs. (58) and (59) 

Xj = pi + x p casd 

yj =P2 + Xpsin 8 

where (pl,p2) is the vector connecting the centers from jth to pth cracks expressed in 
the jth coordinate system: 


Pi = ( r pY ~ r fy) sin (fij + (r p x — r jx) cos Vi 
P 2 = (r p y - rjy ) cos (pj - (rpx - r jX ) sin <pj 


(62) 


The regular, normalized form of the parameters Q f for * — 0,1,... ,6 is obtained by 
using the coordinate normalization x p = Op£ and tj = ajT, in addition to the coordinate 

transformation, to produce the parameters Q» : 


Q< resr) = oc [(a + c) 2 + (6 - d) 2 ] [(a/r - Pl - Opf cos 0) 2 + 

26(a,r - pi - <*,£ cos 0) (p 2 + <Zp£sin0) + (a* + 1?) 0>2 + <*pf sin0) ] 

[(a, T - pi - ap( cos e ? + 2d(°i r -Pi-a^cos0)(p2 + a,^sin0) + 

(c’ + ^fo + Op^) 2 ] 

Qf-es) _ Rl { ajT - Pi - Opt cos of + (P 2 + Op£ sin 0) [R 2 (o,r - Pi - Op£ cos 0)‘ 

+ (P 2 + <*p£ sin0) R 3 (tyT - Pi - ®p£ cos0) 4- (?2 + %£ sin0) 2 i^] 

Q(reg) _ f( 0iT - Pl - Op£ cos0) 3 Re + (jp2 + Op£ sin 8) [r 7 ( ajT - Pi -a^caaO) 
- (P 2 + apt sin 8) Rg (a,r -p x - Op£cos0) + (pa + <Jp£ sin0) 2 f?»] } i?s 


(63) 


(64) 


(65) 




Rio ( aj r -pi-apZcosdf + fa + Opt sin 8) [- (o,r - pi - Op£ cos 0) 2 R n 
+ (P 2 + Op£sin0) (ajT - p\ - apt cosd) Ria + (P 2 + a,4 sin0) R®] 


( 66 ) 


Q ( r ] = 




{(a jT - pi - cos 0) 3 Hu + (P 2 + Sp^sin 0) [Ris (<*;T - Pi “ Op^cas 8) 
+ (p i + a p £sm8) for - pi - flp£ cos 8)Ri 6 + (pa + Op£ sin 8f Rn ] } Rs 

Ri & { aj T -Pi -Oj£ cos 8) 3 + (P 2 + Opfsinfl) [Rit(ajT - pi -Op£cos 0) 2 
- (P 2 + OptsinO) (ajT - pi - Op£ cos 0)Rw + (pa + Spf sin 0) -Raj 


(67) 


( 68 ) 


Of"* - 


|(a,r — pi — Op^ cos 0) 3 R 2 \ + (P 2 + °p£ sin 0) 
+ (P 2 + flp^sin 0) (o,r - Pi - Op£ cos0) R 23 + 


R22 (ajT - pi - flp^ cos 0) 2 
(pj + <ip£ sin 0) 2 R13] } -^5 


(69) 
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So the regular normalized component of the shear stress becomes 
& 3 XpVp = ^ ken fjiitJdTj + ker 2 


keri _ °i 1 - [_ (Q< r< ■») - Q< res) ) sin 0 cos 0 + (cos 2 0 - sin 2 #)] (71) 

27T buQo^ 

ket2 = a i 1 [_ (Q< re »> _ Q< res) ) sin 0 cos 0 + <?S res) (cos 2 0 - sin 2 *)] ( 72 ) 

27r bzzQcT*' 

The regular normalized component of normal stress is 


j 

C ker 3 fjx{tj)dTj + f ker 4 fjiftfidTj 
-1 J ~ l 

( 73 ) 

a i 1 

sin* 9 + COS 1 9 - 2Qi”* ) sin 9 cos 9] 

( 74 ) 

^bnQ^' 


aj 1 

[Q? M) sin 2 9 + Qi" 9 ’ cos 3 9 - 20^”’ sin 9 cos 9] 

. « r 

( 75 ) 

27T 

>1 * — 


uuwj \ yy • 

forau ^ _ f ji^ x f n dT + Jij ker 2/12^ + ... + J- 1 ker 1 /(n-i)i rfr (76) 

+ fix ker 2 /(n-i)2dr + ^ fix fydr + *■- /-1 £| dr } 

ker 3/nrfr + fix ker 4 / 12 dr + ... + /- 1 ker 3 /( n -i)idr (77) 

. , , ri ikL dr + ^ 

+ /_ x ker 4/(n^l)2aT + -7- J-l r-e° T ^ » J “ 1 r ~* J 

The formulation of this system of singular integral eqmrtions is complete > ««** 
single -value conditions for the auxiliary functions /„ are chosen. In the case of str g 
cracks, this single-value condition [ 3 ] is 

j'jin(r)iT = 0 < 78 > 

where j stands for the jth crack, and tj takes on the value of 1 or 2. 

5 SOLUTION FOR THE STRESS INTENSITY FAC- 
TORS 

The integral equations obtained are of the Cauchy type; thrn for sharp cr^ fe stresses 

and strains Z have a square-root singularity, «d ^ “ 
used (see refs. [ 4 ], [ 5 ], (6) and [ 7 ]). Therefore, the mode I and II SIF s for the yth 
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JtJ(l)= Mm|2(£ - 1)]* {,<£(£, 0)} (79) 

*5(1) = lim[2(£ - l)]i {jojK.O)} (80) 

*5(-l) = Bm[-2(1 +£))i {,<£«, 0)} (81) 

*5(-l) = f Um i [-2(l + £# {,«£,(£,<>)} (82) 


where the nomial end shear stresses, eqs. (76) and (77), are used. Note that this 
definition of a SIF was originally developed for an isotropic material wherein the mode- 
I normal stress is related only to the normal crack-opening displacement At) and the 
mode-II shear stress is related only to the shear displacement Au. The same definition 
can also be applied to anisotropic materials by assuming that modes I and II are based 
on normal and shear stresses only; however, as will be shown in eqs. (86) to (89), these 
stresses are driven by a mixed mode displacement (Av and Au) field. 

It is well known [3] that the auxiliary functions (/) can be expressed as a product of 
the unknown bounded functions (H) and the known singular weight functions (w): 


/(r) = H(t)w{t) 

The singular weight function w for a sharp crack is 

(83) 

«»(r) = (r 2 - I)"* 

(84) 


Erdogan [3] found, for example, that in the case of a Cauchy-type singular integral 
equation (eqs. (76) and (77)), the dominant part can be expressed in terms of the 
function H evaluated at the tips of the jith crack: 


w J\ = H,i(-1) $3r (T + 1)_i _ ^ (1) ^ (T - ir * + 0<T) (85) 

where rj is 1 or 2, and 0(r) is the higher order term, which in subsequent calculations 
is neglected. Equations (85) can be substituted for the dominant part (last two terms in 
eqs. (76) and (77)) of the normal and shear components of the total stresses in eqs. (79) 
to (82). This substitution and subsequent evaluation of the limits at the crack tips results 
in the redefining of the SIF’s, normalized with respect to y/a[ and a° jp , and expressed in 


terms of the functions H n j". 

JfeJ(l) = [E 3 i l) Hij( 1) + E* 2) H V ( 1)] (86) 

H(l) = [j^tfiiCl) + Bg>H v { 1)] yfcfa ( 87 > 

jbi(-i) = - [e’A 1 ) h v {-i) + m 
kU- 1) = - [E^Hvi-l) + EffH v (-1)] sfchi (89) 


The Lobatto-Chebyshev collocation integration technique was employed, because it 
is known to provide excellent results in dealing with the preceding Cauchy-type singular 
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intsgr.1 equations. The unknown functions H v , are determined at a dis^ete set of points 
n r 2 r m called abscissas. In this way, each integral equation is reduced to a set o 

afeetaaic equations with unknowns ff„(n), ftn(r,) ftu(r m ), whiA are Uie Ascmte 

values of the functions hence its name, discrete auxiliary function. Notethat H Xj an 

H 2j are proportional, respectively, to the difference in shear and normal displacements 

at the crack tips: 

Au ~ Hi (90) 

Av ~ Hi 

Consequently, tf ni can be used as a measure of the crack-opening displacement at the 

“tSTof the singular integral equations subjected to the stress boundary condition, 
(eqs. 5) can be replaced by m - 1 algebraic equations with 2nm untaown parameter, 
(see ref. [2]).In the Lobatto-Chebyshev method, the abscissas are calculated according 

to 


r r = cos — — ^ for r = 
m — 1 


with the corresponding weights given by 

7T 


Wl = W m = 


and w r = 


7T 


for r = 2,3,...,m- 1 


(91) 


(92) 


2(m— T) m — 1 

The collocation points are then found by using the formula 

£. = cos — — for z — 1,2, 1 

2m — 2 

In order to have the complete system of 2nm algebraic equations, the single-value con- 
ditions (eqs. (78)) are also expressed by using the collocation technique: 


(93) 


m 

£)«u(T r ) UV = 0 


(94) 


r= 1 


Thus, the resulting system of algebraic equations can be written in the form 

[A]{H} = {ft} ( 95) 

where [A] is a fully populated 2nm x 2 nm matrix of coefficients and {11} is the loading 

^ThTuS^o^n parameter vector {H} can be determined through inversion of the [A] 
matrix; thus, 

(Jf) = [A]-'{K} t 96 * 

although only the appropriate values (<.«., ff_(± 1)) are used to calculate the SIF's for 

the ith crack (see eqs. (86) to (89)). , , .. 

The general solution for any multicrack problem is now complete with the autonmtac 

generation of the associated FORTRAN code for the evaluation of eqs. (96). Tins 
FORTRAN program was utilized to obtained the following results, which are compared 
with published results obtained by other methods. 
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6 STRAIN ENERGY RELEASE RATE 


In fracture mechanics, perhaps the most important physical quantity is the strain energy 
release rate (G r ), usually denoted by G. Cherepanov[8] discussed the generalized formula 
for G t given an anisotropic material. 


G = i JT* {<&(*, 0) {v (x - da, 0 + ) — v(x — rfa.O - )] + 
o£ (x, 0) [u (x - da, 0+) - u (x - da, 0')]} dx 


(97) 


Using the roots of the characteristic equation in terms of the real components eq. (11), 
we can show that 


G t = - 


(a + c)(ac -M) + (b + d)(ad + bc) ^ + fc|(o + c) ^ 
(ac — bd) + (ad 4* 6c) 


(98) 


or utilizing eqs. (102) and (105), developed subsequently, we can obtain a simplified 
expression for Gr, that is 

Gt = ( a + c ) 


where 

K *ff ~ 

and represents an effective SIF. Note the significant coupling between the normal and 
shear stresses and the displacement components in this effective SIF. 

hi the case of an isotropic material we may substitute o = c = 1, 6 = d = 0, 
E St ~ E vv> and 6n = bn = i into eq. (98). Consequently, the well-known fracture 
mechanics relationship is recovered: 


p(i) 

I 


<**■>= l {[*n l + k~t} 


( 100 ) 


7 NUMERICAL APPLICATIONS 

The focus of the following parametric study will be limited to investigating the influence 
of the crack geometry configuration and strength of anisotropy on the resulting driving 
force. To accomplish this the four independent elastic constants were taken to be 


Glt 

Err 

vlt 

Ett 

Ett 


0.4 

0.25 

1.0 


and the strength of anisotropy, §££,is a specified constant greater than or equal to one. 
Note that Err is always chosen to be the weaker direction, i.e., to be less than or equal 
to Ell- Although, the influence of Glt and Vlt on the driving force is important, this 
aspect of the parametric study will be reserved for future work. 
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7.1 Two Crack Interaction 


7.1.1 Collinear Cracks 

In order to validate the results obtained with the automatically generated TOKTRAN 
code, the well-known problem of two parallel interacting cracks is consjdered here^The 
plate, with two cracks of length 2a and a preferred direction defined by ’ “ 

subjected to a normal far-field stress state (apy) as shown in Fig. 3. Results are obtained 



Figure 3: Geometry and loading condition of two collinear cracks problem. 

over a wide range of strengths of anisotropy, as defined by the ratio E LL fEjr.AsJ\g. 4 
shows, the SIF’s do not depend on the strength of anisotropy, even though the discrete 
auxiliary functions do. The results indicate that at both the inner and outer crack tips 
mode-I SIF’s are exactly the same as the isotropic SIF’s from references [9] and 111], an 

mode-II SIF’s are zero for this configuration. . 

Although, the mode-II SIF is zero over the entire range of strengths of anisotropy 
examined, the shear crack opening, as represented by H x (see eq. (90) is zero only °r 
the special case of an isotropic material (i.e., Ell/Ett = !)• Figure 4 dearly shows that 
even small amounts of anisotropy (. Ell/Ett > 1) produce shear displacements at the 
crack tip and that this shear displacement increases significantly for 1 < Ell/ Ett 
and becomes constant for E LL /Err > 15. Consequently, two collinear cracks within a 
transversely isotropic material will always (provided ft = 22°) produce a mode-I crack tip 
local stress field with a mixed-mode local displacement field even when the strengtho 
anisotropy is small. This fact can be understood, even for the case of a single cracky], 
by examining eqs. (86) to (89), where it is apparent that in an anisotropic maten 
the normal and shear stresses are coupled with both normal and shear displacement 
components, respectively. 
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Strength of anisotropy, E^/E^ 


qT 
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Figure 4: Discrete auxiliary functions and normalized SIF’s versus strength of anisotropy 
Ell! E rr for two equal collinear cracks in a composite plate reinforced at D = 22° and 
subjected to far-field normal stress, see Fig. 3. 



Figure 5: Modified stiffness parameters versus strength of anisotropy Ell! E rr f° r two 
equal collinear cracks in a composite plate reinforced at fi = 22° and subjected to far-field 
normal stress, see Fig. 3. 
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o 5 10 15 20 25 30 35 40 

Strength of anisotropy, 

Figure 6: Normalized Gr versus strength of anisotropy Eu/Brr for two equal coffinear 
cracks in a composite plate (reinforced at « = 22"; and subjected to far-held nonrnd 

stress, see Fig. 3). 

Figure 5 shows how the MSP’s vary with respect to the strength of anisotropy 
Ell/Ett- Using the isotropic values (i.e., E ll /Ett = 1) of the MSP’s, we can o 
tain the results given in equation (57). Increasing the strength of amwtropy . 
increases 1 B (1) , E& and Eg, while Eg remains nearly constant. Note that E^/ ^yv 
over the entire range of Ell/Ett- Revisiting eq. (87) and Figs 4 and 5, ©plains w y 
the mode II SIF is zero. It stems from the fact that although both normal and shear dis- 
placements are induced, Eg is so much less than Eg that the influences of the normal 

and shear displacements are counteracted. , . . 

As discussed previously, the total Gr represents an important measure of the driving 
force for crack propagation in fracture mechanics. The G T ’s at the inner and outer cr 
tips for the two collinear cracks are shown in Fig. 6. Clearly, the maximum Gr s (for 
both the inner and outer crack tips) occur in the isotropic case. In the anisotropic case 
when E ll /Ett > 15 the G T ’s rapidly reduce to nearly 50% of the isotropic values. 
Although both the Gr (Fig. 6) and the SIF’s (Fig. 4) indicate that the urner crack tap 
will propagate, only the G r ’s unambiguously indicate that as the strength of aniso ropy 
increases, the crack-driving force is reduced. Thus, we can conclude that (1) the isotropic 
case gives rise to the greatest driving force and (2) that the SIF (unlike the G r ) is unable 
to detect the decrease in the crack-driving force as a function of strength of ani«rtropy. 

Now let us examine the influence of changing the preferred direction angle Si on the 
same collinear crack configuration with a strength of anisotropy ratio Ell/Ett - 4U. 
Results shown in Fig. 7 indicate that the mode-I SIF’s are the same as those previously 
obtained, even though the discrete auxiliary functions show that the normal and shear 
crack-opening displacements vary significantly with preferred direction. Note that the 
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Figure 7: Discrete auxiliary functions and normalized SIF’s versus O for two equal 
collinear cracks in a composite plate (strength of anisotropy En,f Ett — 40; subjected 
to far-field normal stress, see Fig. 3). 
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shear displacements are zero only for the special cases in which the preferred direction is 
parallel (ft = 0°) or normal (ft = 90°) to the crack configuration; otherwise both openings 
haw nonzero values. The local shear crack-opening displacements maxima occur at 
ft = 22° and 68°, whereas the local nonzero minimum for the shear displacement occurs 
at ft = 45°. Conversely, the normal crack-opening displacements are apprcoama y 
constant over the range 0 < ft < 40 with the maximum value occurring at ft - dlT; it 
then significantly decreases when ft > 40, until a value equivalent to 30% of the maxmrum 
normal opening is reached at ft = 90°. Clearly, this indicates that although mixed mode- 
displacements may be observed, in actuality, counter to common expectations resulting 
from isotropic materials, the anisotropic case may have only a mode-I driving force 
In Figure 8 the MSP’s are shown as a function of the preferred direction orientation 
angle ft. Note that the E^} curve relative to the curve E$ is symmetrical abou e 
line ft = 45°. The other two MSP’s are again equal to each other and to zero for the 
orthogonal cases ft = 0° and 90°. Consequently, E<$ and must be related to bis 
and has (see eq. (9)). In order to give a clearer physical interpretation to and 
let us relate them to the corresponding material parameters Ell end Ett- 

Considering the fact that . * 

E < 2 > = E^ l) l lul l 


we can substitute eqs. (54) and (55) into eq. (101) and solve for ft** 

&22 = &11 ( fl2 + k 2 ) (c 2 + d 2 ) 

Or in terms of the roots of the characteristic equation 




( 102 ) 


(103) 


In the case of a transversely isotropic or orthotropic material, 6 = d - 0, so eq. (103) 
reduces to a well-known mathematical relation of the following form. 

22 = Ell (104) 

Ett 

By substituting eq. (102) into eq. (56), and the resulting equation into eq. (53), we 

fc “ Uhl “ ^ c(a» + (105) 

~ n (a + c) 

Thus, by using eq. (104) together with the transversely isotropic or orthotropic form of 
eq. (105) (i.e., 6 = d, = 0) we can show that 


= ac — 


b=d = 0 



(106) 


On the basis of eq. (106) and our observations from Fig. 8 we may conclude that 
for the general anisotropic case, is proportional to the square root of the effec ive 
Young’s modulus in the local x-direction and E$ is proportional to the square root of 
the effective Young’s modulus in the local y- direction. Note, however, for the general 
anisotropic case the constants of proportionality are interdependent. 
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Figure 8: Modified stiffness parameters versus (2 for two equal collinear cracks in a 
composite plate (strength of anisotropy Ell/ Ett ~ 40; subjected to far-field normal 
stress, see Fig. 3). 



Figure 9: Normalized Gr versus ft for two equal collinear cracks in a composite plate 
(strength of anisotropy E^jExr ~ 40; subjected to far-field normal stress, see Fig. 3). 
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The importance of being able to use the G r as a damage propagation criterion is 
once again illustrated in Fig. 9, which shows the total G r versus the orientation angle SI. 
Clearly, there are a number of angles at which G T extrema exist. The global maximum 
occurs at Q - 0°, and the global minimum at ft = 90° , for both inner and outer crack 
tips. The other two important angles at which the G r ’s reach a local minimum and 
maxmami are Q = 12° and ft = 42°, respectively. Again a comparison of Figs. 7 and 
9 shows that the Gr’s can indicate the critical angles of interest, whereas the SlF s are 
completely insensitive to variation of the preferred direction. 


7.1.2 Parallel cracks 



Figure 10: Geometry and loading condition defining two parallel crack problem. 


Now lot us consider the variation of the SIF’s, the crack-opening displacements, end 
the Gr’s for the case when the two cracks are not collinear. A convenient parame er 
that can be controUed is the inclination angle designated a between the horizontal axis 
and the line connecting the inner crack tips in Fig. 10. An example of the previous case 
(collinear cracks) is easily obtained by setting a = 0“ . The constant parameters are the 
distance between inner crack tips, O.lo! (where 2a, is the crack length of crack nmnber 
1), the strength of anisotropy Ert/Bir = 40; and the preferred direction Cl = 45. As 
a result of changing the crack configuration, the SIF’s are no longer constant for the 
inner crack tips; those associated with the outer crack tips, however, are not significant y 
influenced. As illustrated in F.g. 11, mode-I SlF for the inner crack tip ' J 

increases with the angle a; whereas mode-II SIF is zero only for a an 


attains a maximum at a = ft — 45°, as we might expect. . , 

It is interesting to analyze the discrete auxiliary functions which approximate the 
shear and normal crack-opening displacements, A u and Av. The shear term of e inner 
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o Shear displacement Au 

• Normal displacement Av Open symboles denote outer tip 
A Mode-1 SIF, Solid symbols denote inner tip 
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Figure 11: Discrete auxiliary functions and normalized SIF’s versus a for two equal non- 
collinear cracks in a composite plate ( strength of anisotropy Ell/ Ett — 40; reinforced 
at ft = 45°; subjected to far-field normal stress, see Fig. 10). 
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crack tip is zero only when a = 40° (even though the shear stress is at a maximum). It 
is larger for this configuration at a = 90° (when the shear stress is zero) than it is for 
the crack configuration. Normal crack-opening displacements are the largest 

also for a = 90° . Note that changing the crack configuration by changing a does not 
change the material properties because they are defined in the local coordinate system; 
therefore, given the preferred direction ft = 45° we can find the pertinent MSP’s in Fig. 

Once A g«»Ti the combined effect of mode- 1 and -II SIF’s and the normal and shear 
local displacements are captured by the single scalar measure known as the G T (see eq. 
(97)). Figure 12 shows the Gr at both the inner and outer crack tips as a function of 
the crack configuration angle a. The shapes of these curves resemble the shapes of ki 
and k 2 for the inner crack tips (the degree of influence that k x and k 2 have on the Gr is 
clearly shown in eq. (99)). Since the maximum values of ki and k 2 are reached at 90° 
and 45°, respectively, the maximum value of the Gt would be expect to occur somewhere 
in between. Figure 12 indicates that for the inner crack tips this maximum is reached 
when a = 70°; however, at the outer crack tip the Gr is only slightly affected (as are ki 
and k 2 ) by a change in a, and it reaches its maximum when a = 90° . 

7.1.3 Nonparallel Cracks 

Consider a related case wherein the material strength of anisotropy is E ll /Ett = 40, 
the preferred direction is denoted by ft = 30°, and the two cracks of equal length are 
configured as shown in the Fig. 13. This figure indicates that crack ab remains horizontal, 
while crack cd rotates around the crack tip c; where the crack tip distance between tips 
b and c remains constant. Here we will examine the influence of varying the angular 
orientation of crack cd, <f>{ 2), from 0 to 180°. Figures 14 and 15 show the mode-I and -II 
SIF’s, respectively, at the four crack tips. Notice that the outer crack tip of the horizontal 
crack (tip a) is only slightly affected by the change in angular orientation of the crack 
cd Whereas, the other crack tips (6, c, and d) display significant and complex interactive 
behavior for both mode-I and -II SIF. For example, mode-I SIF for crack tip b has a local 
mn-yimnm a t 21, 42 and 132°, and a local minimum at 0, 31, 122, and 180°, while crack 
tip c starts at the same value as tip 6, then smoothly decreases to approximately zero 
(within the range 80 < <f>(2) < 140°) whereupon the SIF sharply rises to again the same 
value of SIF as that of crack tip o, at <f>{ 2) = 180°. Similarly, extremum are observed for 
the mode-II SIF at a variety of angle, that is 20, 43, 110, 120 and 132°. 

The actual calculated G r ’s are shown in Fig. 16. By comparing Fig.16 to Figs. 14 
and 15, we can observe that the location of the extremum for the Gt s are similar to the 
SIF’s. However, it is evident that the Gt for outer crack tips a and d are dominated by 
mode-I SIF, whereas inner crack tips 6 and c are significantly influenced by both modes. 
Figure 16 also shows that G r is the largest at crack tip b for 0° < <f>{ 2) < 100° and 
123° < <f>( 2) < 153° indicating possible self similar crack propagation within this range. 
Note that when crack cd becomes aligned with the preferred fiber direction, Gt (at crack 
tip b) reaches a minimum while crack tip c reaches a maximum. Thus, we may conclude 
that at this critical configuration, 4>{2) = 30°, the inner crack tips 6 and c may be driven 
towards each other and connect to form a macro kinked crack. 
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Figure 12: Normalized Gr versus a for two equal noncollinear cracks in a composite 
plate (strength of anisotropy Ell/^tt = 40; reinforced at fl = 45°; subjected to far-field 
normal stress, see Fig. 10). 



Figure 13: Geometry and loading condition problem of two non-parallel cracks. 
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Figure 14: Normalized mode-I SIF’s versus crack-2 angle ^(2) for two equal cradcs m a 
composite plate (strength of anisotropy Ell/Ett = 40; reinforced at Q = 30°; subjected 
to far-field normal stress, see Fig. 13). 



Figure 15: Normalized mode-II SIF’s versus crack-2 angle ^(2) for two equal cracks in a 
composite plate (strength of anisotropy E ll /Ett = 40; reinforced at fi - 30 subjected 
to far-field normal stress, see Fig. 13). 
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Figure 16: Normalized Gr versus crack-2 angle 0(2) for two equal cracks in a composite 
plate (strength of anisotropy Ell/ Ett — 40; reinforced at f2 = 30° ; subjected to far-field 
normal stress, see Fig 13. 



Figure 17: Geometry and loading condition of two inclined cracks. 
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Figure 18: Normalized Gr versus normalized crack tip distance for inclined^cracks in a 
composite plate (strength of anisotropy Ell/Ett ~ 40; reinforced at ft — 30 ; subjected 
to far-field normal stress, see Fig. 17). 

7.1.4 Influence of Crack Spacing 

To determine the influence of crack spacing, let us consider the preceding two crack system 
for the critical configuration $2) = ft = 30° and analyze the behavior of these cracks to 
examine how the Gr varies with respect to the normalized inner crack tip distance d (see 
Fig. 17). As we would expect, given the results shown in Fig. 13, the Gr’s of the outer 
crack tips are hardly influenced by the change in positions of the cracks. But the Gr’s for 
the inner crack tips (6 and c) display a strong interaction, especially for 0 < d < 0.2ai , as 
shown in Fig. 18: the closer the crack tips, the higher the normalized G r , and thus, the 
smaller the far-field stress state required to cause the cracks to propagate toward each 
other. 


7.2 Three Crack Interaction 

Consider a transversely isotropic plate with a strength of anisotropy E ll /Ett = 40 and 
a three-parallel-crack system as shown in Fig. 19. The two cracks denoted ef and 
are always symmetric, with respect to the horizontal line that coincides with crack ah. 
First let us keep the distance between inner crack tips b , c, and e (D/,) constant at 
O.lox, while the angle ft , describing the preferred direction, is varied. Mode-I SIF’sfor 
tips 6, c, and e are shown in Fig. 20. A number of observations can be made from 
Fig. 20. First, the k\ for crack tip 6 resembles an inverted parabola with a maximum 
at ft = 45°. Second, ki for crack tip b is the largest of the three inner tips - for all 
preferred directions because of the magnification influence of cracks ef and cd (situated 
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Figure 19: Geometry and loading condition defining three parallel cracks problem. 



Figure 20: Mode-I normalized SIF’s versus Q for three parallel cracks in a composite 
plate (strength of anisotropy Eli/ Ett — subjected to far-field normal stress, see 
Fig. 19). 
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Fieure 21: Mode-I normalized SIF’s versus strength of snisotropy WErr for three 
parallel cracks in a composite plate (reinforced at fl = 60*; subjected to far-field normal 

stress, see Fig. 19). 

in front of crack ab). Similarly, because of the mutual shi ) ^ 1 ^ *** cracla above 

and below, the k x for cracks cd and ef must be smaller . Thud, the *1 at inner tips 
c and e are identical in the range 0 < ft < 25, but they begin to deviate from one 
another, in a symmetrical manner, for preferred directions ft > 25 Afeo since kj fo 
crack tip e is largTthan that for crack tip c, this would suggest that ^^ence of the 
singular stresses is transmitted over greater distances along the preferred diction, thus 
confirming the concept of stress-channeling (or in our case damage channeling) al g 
the preferred direction, as discussed by Spencer[12]. The distance of influence is c e y 
dependent on the strength of anisotropy specified, as can be seen m Fig. . 

Mode-II SIF’s for the three inner crack tips are shown m Fig. 22. Here, the a 
value of fc 2 for crack tip b is approximately zero for the various preferred directions an 
exactly zero for the orthogonal conditions (0 = 0, 90°) . The absolute v^ueof ^ fo 
crack tip e is the largest of the three with a local maximum at ft = 25 ■ Again, this is 
a function of the directional stress channeling effect. Also, note that the absolute ues 
of k> for crack tips c and e are identical when the cracks are paraM and nonn^ to toe 
preferred directions; however the magnitude at ft = 0 is more Hum twice that for ft - 90. 

The variation of the G T as a function of preferred direction, for toe three mner crack 
tips 6, c, and e are shown in Fig. 23. Clearly, the G T combines all of toe af ^ rem ® n ^^ 
characteristics for each crack tip into one convenient parameter that exhibits * strg 

dependence on the preferred direction angle ft. Local ma ^f^ th ^ CUrV ^ 
within a 10° range centered at ft = 45°. It may be concluded from Fig 23 that crack 
propagation is easiest when the preferred direction makes a 45° angle, th ^ c ^^ mg 
through reinforcement (or damage channeling) crack tips b and e. Hence, as a result of 
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Figure 22: Mode-II normalized SIF’s versus angle £2 for three parallel cracks in a com- 
posite plate (strength of anisotropy Ell/ Ett ~ subjected to far-field normal stress, 

see Fig. 19). 



Figure 23: Normalized G T versus £2 for three parallel cracks in a composite plate (strength 
of anisotropy Ell/ E rr — 40; subjected to far-field normal stress, see Fig. 19). 
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Figure 24: Normalized G r versus horizontal component of the tap distance D fc for 
parallel cracks in a composite plate (strength of anisotropy E ll /Ett = 40; reinforced at 
£2 = 45°; subjected to far-field normal stress, see Fig. 19). 

a physical crack propagation, we would anticipate a zigzag crack that develops through 

the connection of cracks ab and ef . .... c » 

Last, let us consider the case of two parallel cracks ef and cd that are a fixed vertical 

distance apart (D. = 0.1a,) and crack ab is slid in between them. Figure 24 ehowstb. 
variation of the Gr’s for all crack tips as a function of horizontal position £>* for the case 
when fi = 45°. When the parameter D h is zero, tip 6 is on the same vertical line tha 
connects tips c and e. Thus, when crack ab is away from the parallel cracks D h > , an 
when tip b is between cracks cd and ef , D h < 0. Note that as cradc tip b comes closer 
to the vertical line connecting tips c and e, the G r s are magnified; with the amphfica ion 
factor of the inner crack tips being significantly greater than that of the outer crack tips. 
Conversely, when D h < 0 , all inner tips become strongly shielded, so the inner Gr s 
sharply drop off, almost to zero. Note that crack tip e has the most rapid decrease and 
quickly reaches a value less than that of tip c. 

7.3 Horizontal Notch Interaction With Three Microcracks 

The final problem to be consider in this paper consist of a transversely isotropic plate, 

with a strength of anisotropy {Ell! Err) ec l ua ^ to 4 ^, OI ^° n _ 

notch (of length 2a x ) and three radially oriented microcracks (of length Oj - a 3 -a A 
O.laj), as shown in Fig. 25. The radial distance between the inner tips of the notch “id 
the three microcracks remains fixed at 0.005a!, while the orientation (angle £2) of the 
preferred direction is varied from 0 to 180°. The Gr’s for all the inner microcrack and 
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Figure 25: Geometry and loading condition defined for a horizontal notch and three 
microcrack problem. 

notch tips are shown in Fig. 26. 

Evidently, the G T for notch tip b is larger than that for any of the three microcracks 
except when the orientation is within two small regions: 8° < ft < 20° and 70° < ft < 98 . 
At these values the Gr of the inner tip of the 45° microcrack ef becomes greater than that 
of the notch. All Gr curves display strong and complex dependence cm the orientation. 
Numerous local extrema exist with the maxima being noted as follows: for notch tip b, the 
local maxima occur at 0, 38, and 140°; for microcrack tip g, a maximum occurs at 30°; for 
microcr ack tip e, maxima occur at 15, 80, and 155°; and for the horizontal microcrack tip 
c, the maxima occur at 0°, 40°, and 150°. Again, if the maximum-Gr criterion is used to 
predict crack propagation, we may conclude that the crack will kink by connecting with 
crack ef for the preferred direction when 8° < ft < 20° and 70° < ft < 98°; otherwise it 
will propagate in a self-similar manner. 

Now let us orient the preferred direction ft at a fixed angle of 15° and vary the radial 
inner tip distance between notch tip 6 and the three microcrack tips c, e, and g. Figure 
27 shows the amp lification effect (which becomes noticeable when the tip distance is less 
then O.loi) resulti ng from the complex interaction of the cloud of microcracks with the 
larger notch crack. Clearly, the G T for notch tip 6 is the largest until the tip distance 
is decreased to approximately 0.04oi, whereupon the G t for the microcrack inclined at 
45° drastically increases - exceeding all other curves and creating the condition for the 

crack to kink. > 

Finally, let us discuss the influence of the strength of anisotropy, as shown in Fig. 28. 
Here, we assume that the preferred direction and radial tip distance are held constant at 
ft = 15° and 0.005ai, respectively, while the strength of anisotropy Ell! E rr is varied - 
When the material is isotropic (i.e., Ell! Ett — 1) tfr e Gr at notch tip b far surpasses 
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Figure 26: Normalized G r versus ft for notch and three microcrack inner tips in a 
composite plate (strength of anisotropy E ll /Ett — 40; distance, 0.005a! ; subjec 
to far-held normal stress, see Fig. 25). 



0 0.05 0.1 0.15 0 2 0.25 0 3 0.35 0.4 0.45 0.5 

Normalized tip distance 


Figure 27: Normalized Gr versus normalized tip distance between inner tips of notch and 
three microcrack in a composite plate (strength of anisotropy Ell/ E hr 40; rei ore 
at ft = 15°; subjected to far-field normal stress, see Fig. 25). 
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Figure 28: Normalized Gy versus strength of anisotropy Ell/Ett f° r inner tips of notch 
and three microcracks in a composite plate (reinforced at f2 = 15°, tip distance, 0.005oi| 
subjected to far-field normal stress, see Fig. 25). 

the Gy’s of the r emaining micro crack inner tips. However, even a small change in the 
strength of anisotropy significantly influences the value of the Gy for all inner crack tips. 
Note that when Ell/Ett > 18, the Gy for the inclined microcrack tip e becomes the 
largest. And since the crack having the highest Gy will propagate first, a kinked crack 
will be generated by connection of crack ef to notch ab. 


8 CONCLUDING REMARKS 

A rigorous formulation has been presented for calculating the crack-opening displace- 
ment, SIF and Gr at the various crack tips of a multicracked anisotropic medium. This 
formulation has been shown to simplify exactly to our previous isotropic formulation 
which was validated for a number of published crack orientations. The size, orientation, 
and distribution of all cracks were considered to be independent parameters of the solu- 
tion. This unique formulations is computationally efficient and offers accurate solution 
capability. It allows us to easily perform numerous parametric studies to analyze the 
contribution of each parameter on the local stress field and the characteristics of the 
damage progression in an anisotropic (e.g. , transversely isotropic) material. 

The problems of two and three, collinear and non-collinear, interacting cracks were 
examined. By varying the strengths of anisotropy, we showed that materials with pre- 
ferred off-axis directions relative to the applied load produced highly mixed-mode crack 
propagation, even when only mode-I type crack geometry was present. A small change 
in the strength of anisotropy (when 0 < Ell/Ett < 5) was shown to highly influence 
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the Gx and crack-opening displacements for a given tip. These parameters then reached 
a plateau when Ell] E rr > 15. Consequently, even slightly anisotropic materials should 
be analyzed by using the fully anisotropic approach discussed in this paper. 

Modified stiffness parameters (MSP’s) were presented as a function of strength of 
anisotropy and preferred direction. We concluded that Ej$ and E^j are proportional 
to the square root of the effective Young’s module in the local x— and j/-direction, 
respectively, and that the other two MSP’s are identical (i.e., Ej$ = E$) and related to 
big and bx such that they vanish under orthotropic and isotropic conditions. 

The discrete auxiliary functions were shown to be related to the crack-opening dis- 
placements A u and Av. We showed that for cases in which the preferred direction and 
the applied stress do not coincide (i.e., off-axis orientations), a mode-I local stress field 
is produced under mixed-mode local deformations, or alternatively, mode-I normal de- 
formation results in a mixed-mode local stress field. Furthermore , the total Gy was 
shown to be the most complete anisotropic fracture parameter (because of its sensitivity 
in detec ting changes in the strength of anisotropy and preferred direction). Thus we sug- 
gest that for an anisotropic material it should be used as the crack propagation criterion 
instead of the SIF. 

Interaction effects were demonstrated for all fracture parameters. Amplification of 
the Gr and/or the SIF was shown to occur when cracks were located in front of the 
main crack. Conversely, reductions in the Gy or SIF were observed when shielding of 
a crack, by other cracks located above and/or below it, was present. Finally, stress or 
damage channeling was discovered to play a significant role in the mechanisms that govern 
crack interaction, in that, stresses were channeled along the preferred direction, causing 
nonsymmetric interaction, even in the presence of symmetric crack configurations. 
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